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Ferroelectric properties in magnetic fields of varying magnitude and direction have been inves- 
tigated for a triangular-lattice helimagnet CuFei-2,Gai02 (a:=0.035). The magnetoelectric phase 
diagrams were deduced for magnetic fields along [001], [110], and [110] direction, and the in-plane 
magnetic field was found to induce the rearrangement of six possible multiferroic domains. Upon 
every 60°-rotation of in-plane magnetic field around the c-axis, unique 120°-flop of electric polar- 
ization occurs as a result of the switch of helical magnetic g-vector. The chirality of spin helix is 
always conserved upon the g-flop. The possible origin is discussed in the light of the stable structure 
of multiferroic domain wall. 

PACS numbers: 77.80.Fm, 75.60.-d, 



Magnetoelectric (ME) effect, i.e. magnetic (electric) 
induction of electric polarization P (magnetization M), 
has long been investigated [ij. The most promising ap- 
proach to achieve gigantic ME effect is the magnetic 
and/or ferroelectric domain control in multiferroics, ma- 
terials with both magnetic and ferroelectric (FE) order. 
In general, the formation of domain structure depends on 
the symmetry ,2]; if each domain is characterized by both 
magnetic and dielectric order parameters, ME effects 
are obtained from modulation of domain distribution by 
electric field [E) or magnetic field {H). This strategy 
was first demonstrated on a ferroelectric weak ferromag- 
net (FM) NiaBrOial with persistent 180°-reversal of P- 
vector under 90°-rotation of i?[3]. However, this type 
of ME control has seldom been achieved, because of the 
rareness of similar FE/FM compounds *4"] and the strict 
symmetry restrictions for selective domain switching [2] . 

One important breakthrough was achieved by the dis- 
covery of ferroelectricity of magnetic origin[5l]. In some 
frustrated magnets, spiral spin order is found to induce 
finite P, which is now explained in terms of the inverse ef- 
fect of Dzyaloshinskii-Moriya (D-M) interaction Q. This 
model predicts the coupling between the sign of P and 
spin chirality (left-handed or right-handed spin rotation) , 
which was experimentally confirmed by studies of polar- 



To apply multiferroics to the novel high-density non- 
volatile storage device, such as an iJ-controUed ferro- 
electric memory, the persistent switch of P-direction as 
well as increase of the number of switchable meta-stable 



ized neutron scattering!?^ 
TbMnOai and MnW04 



. In some helimagnets like 
10( 1 , a specific direction of static 
H induces the 90°-flop of spin-spiral plane and accord- 
ingly of P-direction. Some hexaferrites show smooth ro- 
tation of P-vector and spin-spiral plane under rotating 
7j[ll|,[l2|- Likewise, the collinear spin order can also be 
coupled with P, due to symmetric exchange striction S]; 
TbMn205 with collinear spin order shows //-induced 
180° reversal of P[l^. In all these compounds, however, 
the P-direction change is not persistent after the removal 
of iJ, since directional change of P originates from the 
//-induced modulation of spin structure. 



states are highly desirable [IJ]. While several approaches 
have realized the //-induced persistent change of P- 
directionjs, 15, I^j, they are all limited to the 180° switch 
between ±P states. Aside from the non-volatile nature, 
//-induced discontinuous switch of P-direction among 
more than two different axes (i.e. directional switch of 
P by other than 90° or 180°) has never been achieved. 
In this study, we extend the concept of domain switch- 
ing to demonstrate the persistent magnetic control of six 
ferroelectric domains with different P-directions for tri- 
angular lattice ferroelectric helimagnet CuFei_a;Gaj;02 
(a;=0.035). The flop of magnetic modulation vector can 
be induced upon every 60°-rotation of in-plane H around 
the c-axis, which leads to every 120°-flop of P-direction 
within the triangular-lattice basal plane. Interestingly, 
the chirality of spin-spiral is always conserved upon the 
P-flop. We discuss the possible origin in the light of the 
stable structure of multiferroic domain wall. 

CuFe02 is known as a member of triangular lattice an- 
tiferromagnets with spin frustration. Magnetic moment 
is carried by Fe'^"'' ion with S = 5/2, and each element 
(Cu/Fe/0) forms triangular lattice respectively. They 
stack along the c-axis in order of Cu+-0^^-Fe^+-0^^- 
Cu"*", and crystalize into the delafossite structure with 
centrosymmetric space group P3m 17| . CuFe02 under- 
goes subsequent magnetic phase transitions (CM-4 NC 

CM-5) under H applied along the c-axis (//[ooi]), while 
keeping the in-plane component of magnetic modulation 
vector {q) parallel to (110) [l^,[l^. Hereafter, we treat q 
as the director (with no -I-/- sign). CM-4 and CM-5 have 
collinear spins along the c-axis, and correspond to ttii 
and tttii spin structures, respectively [lijj. Ferroelectric- 
ity is observed only in the NC phase with proper screw 
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FIG. 1: (color online). Temperature (T) dependence of mag- 
netic susceptibility (x) and [110] component of electric polar- 
ization (P). Magnetic field {H) is applied along (a),(b) [001], 
(c),(d) [110], and (e),(f) [110], respectively. Arrows indicate 
the direction of T-scan. In (c), a magnified profile of x a-t 
O.IT (arbitrarily off-set) is also shown, (g)-(i) Three out of 
six possible multiferroic domains with proper screw magnetic 
structure on triangular lattice. Circled "R" and "L" denote 
the chirahty of spin spiral. Directions of P and magnetic q- 
vector are also indicated, (j)-(m) Distribution of multiferroic 
domain(s) favored under various H and electric field (E). The 
spin chirality corresponding to each P-domain is shown in (j). 



magnetic structure [20], where spiral spin rotates within 
the plane perpendicular to the magnetic a- vector pi). 
The P in NC(FE) appears parallel to g- vector (Fig. 
1(g)), and the spin chirality is confirmed to be reversed 



for the reversal of P-direction (Fig. 1 (h) ) [2^, lil H] . The 
present ferroelectricity with P || (? cannot be explained by 
the simple inverse D-M model, which predicts P _L g[6|. 
Instead, the variation of 7r-bonding between Fe and O 
ions under the influence of spin-orbit interaction has been 
suggested as the origin of ME coupling IS 13, 2^ • The 
critical H to induce the NC(FE) phase can be reduced to 
zero by nonmagnetic impurity doping on Fe'^"'" site. This 
was originally shown for Al-doping but with much re- 
duced P H, 27] due to strong pinning of domain wall 22). 
This problem can be overcome by Ga-doping[2^, and 
here we adopt CuFei_a;Gaa;02 (a;=0.035) as a target mul- 
tiferroic. Without external field, this compound has six 
equivalent ferroelectric domains with P ]] (110) (Fig. 
l(j)) due to high symmetry of underlying triangular lat- 
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FIG. 2: H-T phase diagrams with H parallel to (a) [001], (b) 
[110], and (c) [lIO] direction. Circles, triangles, squares, and 
diamonds are the data points obtained from P-T, P-H , x~T, 
and M-H curves, respectively. All the data were taken from 
the increasing T ox H runs[13]. Ferroelectric (FE) state is 
observed in the shadowed region, (d)-(f) H-dependence of P 
(gray line) parallel to [110] and magnetization (M: black line) 
at 2K, with H applied along (d) [001], (e) [110], and (f) [lIO] 
direction. In all measurements, H was swept from OT to 14T 
and then back to OT, after T was lowered to 2K at OT. The 
arrows indicate the direction of Jf-scan. 



tice. From the P || (7 relationship, they correspond to six 
spiral magnetic domains with three different q \\ (110) 
axes and two spin-chiral degrees of freedom. 

Single crystals of CuFei_a;Gaa;02 (x=0.035) were 
grown by a floating zone method [tT]. They were cut 
into a rectangular shape with the paired (110) and (110) 
surfaces, on which silver paste was painted as electrodes. 
To deduce P, we measured the polarization current with 
constant rates of temperature(T)-sweep (2K/min), H- 
sweep (80 Oe/sec), or iJ-rotation (l°/sec), and inte- 
grated it with time. To enlarge the population of specific 
P-domains, the poling electric field [E — 250kV/m) was 
applied along [110] in the cooling process and removed 
just prior to the measurements of polarization current. 
M was measured with a SQUID magnetometer. 

Since only the properties under zero magnetic field 
have been reported so far for Ga-doped CuFe07(28|. we 
first estabfish H-T phase diagrams for the 2:=0.035 spec- 
imen in i?[ooi], -^[110] and i?[iio] (Figs. 2(a)-(c)). They 
are determined from the measurements of T- and H- 
dependence of M and [110] component of P (P[iio]) (Figs. 
l(a)-(f) and Figs. 2(d)-(f)), by analogy with the case 
for CuFei_:,AU02 (a;=0.02) under ir[ooi] [H, [i^. Here, 
paramagnetic and two different sinusoidally-modulated 
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FIG. 3: (color online), (a)-(c) [110] and (d)-(f) [110] compo- 
nents of P simultaneously measured in H rotating within the 
(001) plane. 6h denotes the angle between //-vector and the 
[110] axis (see Fig. 4(d)). Arrows indicate the direction of 
//-rotation. Absolute value of P was determined by T-scan. 



collinear incommensurate magnetic phases 3 2^ are re- 
ferred to as PM, ICMl and ICM2, respectively. With 
any direction of H, the boundary of the FE phase (shad- 
owed region) always coincides with that of the NC mag- 
netic phase, which ensures the coupling between ferro- 
electricity and proper-screw magnetic structure. In T- 
scan profiles, the onset of spiral magnetic order, coupled 
with the emergence of ferroelectric P, can be detected as 
the sudden drop of magnetic susceptibility x{— M/H). 
While the NC(FE) phase is replaced by CM-5 under 
H[ooi] > 12T, we found that NC(FE) survives against 
in-plane H (//[no] and //[iio]) up to 14T. 

Even with in-plane //, the magnitude of ^'[iio] shows 
significant //-dependence, as shown in Figs. 1(d) and 

(f): -P[iio] increases with //[no] and decreases with //[j^jo] • 
This behavior can be interpreted as the rearrangement 
of multiferroic domains, as described below. In general, 
antiferromagnetically ordered spin moments prefer to lie 
within the plane perpendicular to //, as typically ob- 
served for the spin flop transition. Thus, in case of proper 
screw magnetic structure, H favors the magnetic domain 
with q II H. In contrast, electric field E affects the selec- 
tion of spin chirality 23[. U E is applied along [110], three 
multiferroic domains are selected as depicted in Fig. l(k). 
Further application of H should sort domains; a domain 
distribution as illustrated in Figs. 1(1) or (m) is favored 
with //[no] or //[iio]- Comparing these arrangements 
(Figs. l(k)-(m)), P[iio] should increase with //[no] and 
decrease with //[iio] • This idea well explains the observed 
results, and similar effects have recently been reported 

Mm- 



for multiferroic delafossite CuCr02 
and //[ijo]"Scan at 2K, both M and P[iio] 
lies around 5T only in the field- increasing run (Figs. 2(e) 
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FIG. 4: (color online), (a) Relationship between the direc- 
tions of P and H , both confined within the (001) plane. 9p 
(Oh) denotes the angle between P- (//-) direction and the 
[110] axis. (b),(c) Corresponding variations of M. Dashed 
lines indicate theoretically expected behaviors, (d) Relation- 
ship among P, q and H a,t Oh = (60n)° (n: integer). 



and (f)), suggesting the field- irreversible rearrangement 
of domain distribution at this specific field direction. 

Next, we have investigated the vector components of 
P in response to the in-plane H rotating around the c- 
axis. For this purpose, the both [110] and [lIO] compo- 
nents of P were measured simultaneously with two pairs 
of electrodes. In this configuration, both P and H can 
be expressed as the two-dimensional vector on the (001) 
plane. Hereafter, we define the angle between H (P) and 
the [110]-axis as 9h (dp)- Since the specimen was cooled 
with H and E both apphed along [110], we assume the 
uniform initial domain state as shown in Fig. 1(1). 

Figures 3(a) and (d) show P[iio] and f [iio] as a function 
of 9h, measured at // = 6.5T without E. Both P[iio] 
and -P[iio] show a periodic change with the cycle of 180°. 
To see the development of P more directly, we plot the 
obtained 9p against 6h for the ^//-increasing run (Fig. 
4(a)). In agreement with the expected initial state in 
Fig. 1(1), the relationship P || // || [110] is confirmed at 
6h = 0. As 9h increases, P suddenly flops by about 120° 
at 9h = 41°, and the relation that -P ]] H \\ [100] holds 
at 9h — 60°. Since H favors multiferroic domains with 
q II P II ±//, this transition can be considered as the flop 
of magnetic g- vector from g |] [110] to q \\ [100]. Note that 
both transitions to P ]| [100] and P ]] —[100] states seem 
possible from the P ]] [110] initial state, but in reality 
only the P ]| —[100] state is selected here. Such a g-flop 
(or 120°-flop of P) is observed upon every 60°-rotation of 
H, consistent with the symmetry of underlying triangular 
lattice. The absolute value ]P] (not shown) is almost 
constant for any 9h- 

Between the two opposite directions of //-rotation, a 
relatively large hysteresis is found in P. If we estimate 
the critical 9h from the average of both direction of scans, 
the g-flop is always centered at 9h = (30 -I- 60n)° {n: in- 
teger). This agrees with the equilibrium point of two 
magnetic g-vectors. The appearance of hysteresis means 
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that an excess gain of Zeeman energy is needed to over- 
come the potential barrier height. When H is reduced 
from 6.5T to 4.5T, the reduction of Zeeman energy leads 
to expansion of the hysteresis (Figs. 3(b) and (e)). Be- 
low 3.5T, the potential barrier cannot be overcome and 
no g-flop behavior is observed (Figs. 3(c) and (f)). 

Corresponding variation of M as a function of Oh is 
shown in Figs. 4 (b) and (c). To obtain the q \\ [110] 
initial state, the specimen was cooled at Oh = 0° with 
H =^ 6.5T and then magnitude of H is fixed prior to 
measurements. Without g-flop, this initial state should 
give AM oc (x|| — X±) cos(20//), where X|| and X|| denotes 
X parallel or perpendicular to q. This agrees well with 
Oh dependence of M at 3T (Fig. 4(c)), indicating the 
robustness of single-g state and the non-volatile nature 
of g-domain distribution. In contrast, the profile at 6.5T 
has a period of 60° and corresponding sinusoidal curve 
is shifted by 60° upon every P-fiop transition. This con- 
firms the emergence of g-flop as the origin of P-flop and 
the clamping of ferroelectric and magnetic domain walls. 

Figure 4(d) illustrates the relationship between P and 
H at 6'jy=(60n)°. At each Oh, we could confirm P || iff, 
in agreement with the g-flop model. Importantly, upon 
each transition, P \\ H and P \\ —H alternately ap- 
pears. Generally, magnetic domains can be mutually con- 
verted by symmetry operation that is broken by magnetic 
orderQ. If we apply space inversion to Fig. 1(g), a do- 
main with opposite P and reversed spin chirality can be 
obtained (Fig. 1(h)). Mirror operation on Fig. 1(g) gen- 
erates another domain with reversed spin chirality (Fig. 
l(i)). Likewise, we can reproduce all six P-domains and 
determine their corresponding spin chirality (Fig. l(j)). 
From this relationship, it is concluded that the chirality 
of spin spiral is always conserved upon the g-flop. 

Recently, a similar systematic behavior of spin chiral- 
ity has been observed in some ferroelectric helimagnets: 
Upon 90°-flop of g under rotating H on ZnCr2Se4, the 
spin chirality is preserved in lower- 77 region but reversed 
in higher-ff region ■ 3 2j. In case of 90°-flop of spin spi- 
ral plane (with fixed g) in unidirectional H on MnW04, 
spin chirality after the transition can be selected by slight 
tilt of iJ[33|. Since two opposite chiral states are ener- 
getically degenerated under H, selection of odd chiral- 
ity upon magnetic transition was explained by the en- 
ergy difference between two possible domain wall (DW) 
structures connecting domains with the same or oppo- 
site chirality. The behavior of spin chirality depends 
on the given ff-condition, therefore the stability of DW 
should be modified by H . This situation likely occurs 
in CuFei_a;Gaa;02: When rotating H induces the P-flop 
from Pi to P2, the angle between Pi and P2 across the 
DW becomes 120° for the same chirality (120°-DW), and 
60° for the opposite chirality (60°-DW). The present ro- 
bustness of spin chirality can be explained, only provided 
that 120°-DW is more stable than the 60°-DW. Detailed 
theoretical calculation of relative stability of two types of 



DWs considering all magnetic, dielectric, structural, and 
chiral degree of freedoms, as well as direct observation of 
domain dynamics, is highly desirable. 

In summary, we demonstrate the persistent magnetic 
control of six ferroelectric domains for triangular-lattice 
ferroelectric helimagnet CuFei_a;Gaa;02 (a;=0.035). The 
flop of the magnetic modulation vector is induced by 
every 60°-rotation of in-plane H around the c-axis, 
which leads to every 120°-flop of P-direction within the 
triangular-lattice basal plane. The chirality of spin-spiral 
is always conserved upon the P-flop, which may re- 
flect the stability of the specific multiferroic domain wall 
structure. In more general, the nature of multiferroic DW 
in a ferroelectric helimagnet should play a key role in de- 
termining the P-direction or spin-chirality upon the H- 
induced g-flop; this may enable the unusual selective P- 
domain switching with varying direction and magnitude 
of H . Ferroelectric helimagnets with highly symmetric 
crystal structure like cubic, tetragonal or hexagonal lat- 
tices can host the multiple 7J-switchable P-domains as 
revealed in this study, and hence are promising for ME 
control with non-volatility and multiple-valued nature. 
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